Quantum theory of fermion production after inflation 
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We show that quantum effects dramatically enhance the production of fermions following preheat- 
ing after inflation in the early Universe in the presence of high excitations of bosonic quanta. As a 
consequence fermions rapidly approach a quasistationary distribution with a thermal occupancy in 
the infrared, while the inflaton enters a turbulent scaling regime. The failure of standard semiclas- 
sical descriptions based on the Dirac equation with a homogeneous background field is caused by 
nonperturbatively high boson occupation numbers. During preheating the inflaton occupation num- 
ber increases, thus leading to a dynamical mechanism for the enhanced production of fermions from 
the rescattering of the inflaton quanta. We comment on related phenomena in heavy-ion collisions 
for the production of quark matter fields from highly occupied gauge bosons. 
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The precise understanding of phenomena out of equi- 
librium plays a crucial role for our knowledge about the 
primordial Universe. Important examples are the gener- 
ation of density fluctuations, nucleosynthesis, or baryo- 
genesis, with the latter being responsible for our own 
existence. 

Most dramatic processes such as far-from-equilibrium 
particle production at the end of inflation pose a par- 
ticular theoretical challenge [l[ . The corresponding non- 
perturbative phenomenon of preheating is conventionally 
described using classical approximations for the bosonic 
inflaton 2]. Their validity for macroscopic occupation 
numbers has been verified explicitly in quantum field 
theory Q. Much less is known about fermion dynam- 
ics in the nonperturbative regime of high Bose occupa- 
tion numbers. Since identical fermions cannot occupy the 
same state, their quantum nature is highly relevant and 
a consistent quantum theory of fermion production after 
inflation is of crucial importance. Related questions arise 
in the context of heavy-ion collision experiments, where 
the production of quark matter fields from highly occu- 
pied gauge bosons can occur far from equilibrium 

So far, preheating dynamics with fermions has been 
mainly investigated based on semiclassical descriptions 
using the Dirac equation with coupling to a homogeneous 
inflaton field |5| . Also backreaction of fermions onto in- 
flaton dynamics has been included. A similar semiclas- 
sical prescription has been used to compute quark pro- 
duction from classical gluons in the context of heavy-ion 
collisions Q. Further inclusion of quantum corrections 
is complicated by a secular perturbative time evolution 
which becomes rapidly invalid. Suitable approximations 
may be based on the two-particle irreducible (2PI) effec- 
tive action, which are known to describe the approach to 
thermal equilibrium Q . Simplified discussions with these 
techniques, which evaluate the dynamics for vanishing in- 
flaton field amplitude, exhibit new nonlinear phenomena 
such as instability- induced fermion production [§]. 

In this work we present a quantum theory of fermion 



production following preheating after inflation. We con- 
sistently include quantum corrections to next-to-leading 
order (NLO) in the Yukawa coupling between the infla- 
ton field and massless fermions. Even for weak couplings 
this turns out to change semiclassical or leading-order 
(LO) results so dramatically that we consider a comple- 
mentary nonperturbative method for comparison. It is 
based on lattice simulations following the techniques of 
Ref . Q . This treats the fermions exactly, but the inflaton 
dynamics becomes classical statistical. In 3 + 1 dimen- 
sions this is computationally expensive and becomes fea- 
sible with the implementation of "low-cost" fermion al- 
gorithms fiol ] . Remarkably, higher-order corrections turn 
out to leave the NLO results practically unchanged for 
the considered range of weak couplings. 

In inflationary cosmology, the Universe at early times 
expands quasiexponentially in a vacuumlike state. Dur- 
ing this stage of inflation, all energy is contained in a 
slowly evolving inflaton field. Eventually the inflaton 
has to transfer its energy to particles, thereby starting 
the thermal history of the hot Friedmann universe. If 
the creation of particles is sufficiently slow (for instance, 
if the inflaton is coupled only gravitationally to the mat- 
ter fields) this transfer may be described by perturba- 
tive decay. However, for a wide range of couplings the 
particle production is governed by a nonequilibrium in- 
stability These instabilities are known to lead 
to exponential growth of inflaton occupation numbers 
in long wavelength modes on time scales much shorter 
than the asymptotic thermal equilibration time. This 
is followed by a turbulent phase with different univer- 
sal scaling regimes for nonperturbative long wavelength 
modes and perturbative higher momenta 12]. De- 
spite the differences between early-universe dynamics and 
the "little bang" generated in a heavy- ion collision, there 
are remarkable parallels concerning the role of instabili- 
ties and turbulence for the process of thermalization [3] ■ 

To reveal the production of matter fields from inflaton 
decay requires the inclusion of fermions. We consider a 



2 



quantum field theory with interaction of the Yukawa type 
- jf^i ((1 ~ 7 5 )4 + (! + 7 5 )^) ^ , (1) 

for i,j = 1, . . . , Nf massless Dirac fermions, where sum- 
mation over repeated indices is implied. (We use nat- 
ural units with h = c = k R = 1.) The interaction is 
constructed to couple the left- and right-handed fermion 
field components ip L = (1 — -f 5 )^/2, ip R = (1 + j 5 )ip/2 
for a nonvanishing expectation value 



<t>(t) = (Tr$(t,x)) 



(2) 



of the inflaton field. This describes the generation of an 
effective fermion mass by spontaneous symmetry break- 
ing. Together with the scattering or decay processes im- 
plied by ([T]), the interaction contains the major building 
blocks required to discuss fermion production. 

To be specific we consider two flavors, Nf — 2, with 
symmetry group SU L (2) x SU R (2) ~ 0(N S = 4), where 
$ = (<t + ittt) /2 with Pauli matrices f and an inflaton 
potential V = to 2 (ct 2 +tt 2 )/2+A (er 2 + tt 2 ) 2 /41N S charac- 
teristic for a linear sigma model. For our purposes it will 
be sufficient to consider dynamics without expansion. In 
chaotic inflationary models, after the end of inflation the 
inflaton field coherently oscillates around the minimum 
of its potential with a very large amplitude. The cou- 
pling of the field to its own quantum fluctuations leads 
to the phenomenon of parametric resonance [l[ . For re- 
alistic scenarios the inflaton mass is required to be about 
6 orders of magnitude smaller than the value of the in- 
flaton field, and the inflaton quartic coupling is about 
10~ 12 . Since simulations are hard to perform for realistic 
hierarchies, we will consider numerical results where the 
mass is an order of magnitude smaller than the inflaton 
with a coupling A = 0.1. The realistic case can then be 
inferred from our analytic results for time scales and pro- 
duction rates. Though all results will be for parametric 
resonance in chaotic inflation scenarios, we concentrate 
on characteristic properties which are independent of the 
underlying instability type and details of the model. 

Fermion production can be computed from 

(3) 

Here F(x,y) denotes the commutator or statistical two- 
point function, and p(x, y) is the anticommutator of two 
fermion fields or the spectral function. In thermal equi- 
librium these would be translation invariant and related 
by a fluctuation-dissipation relation in Fourier space: 
F eq (w,p) = — z [1/2 — nFD(w)]p C q(w, p) with the Fermi- 
Dirac distribution 7Jfd(^) = 1 / '[exp(w /T) + 1] for zero 
net fermion charge [7(. The question of fermion pro- 
duction out of equilibrium may be discussed in terms 
of a time-dependent occupation number n^(t, p), which 
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FIG. 1: Time-dependent fermion occupation number. 



starts out at n^{t = 0,p) ~ after inflation and ap- 
proaches a constant distribution in equilibrium. Calling 
M^,(t) = g4>{t)/Nf we consider a standard definition [B[ 



n^(t, p) 



1 \p\F v (t,t;p)+M^t)F s (t,t-p) 



v /p 2 + M|(t) 



(4) 



where F s = tr(F/4), F v = tr(p- jF/A) with the trace 
acting in Dirac space. These are the nonvanishing com- 
ponents of F at initial time consistent with symmetries. 
The standard semiclassical (LO) evolution equation is 



N 



F(x,y) = 0. 



(5) 



Figure [TJ shows the corresponding LO result (dotted line) 
for the fermion occupation number for £ = g 2 / \ = 0.1 at 
fixed time t = 250/ooMo in units of the rescaled initial in- 
flaton field <fio = 4>{t = 0), i.e. co-Mo = (j>o/^6N s /X with 
inflaton mass Mo . Since in this approximation the result 
is strongly oscillatory, in addition a time average (dashed 
line) over a period of about 20/coAfo is displayed. These 
characteristic LO properties of fermion occupation num- 
bers, which vary periodically in time governed by reso- 
nant excitations, are discussed in Ref. [5j. 

In order to see whether the inclusion of quantum cor- 
rections changes these results, we first employ 2PI effec- 
tive action techniques [7|, 1 1 31 ] . The fermion contribution 
in the coupling expansion of the 2PI effective action at 
NLO in g reads diagrammatically: 



(0) 



Here the dashed lines denote self-consistently dressed 
fermion propagators and the solid line represents the in- 
flaton propagator. For both LO and NLO results the 
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FIG. 2: Total number density of fermions. 



FIG. 3: Inflaton field amplitude. 



same truncation for the bosonic sector is used, which in- 
cludes all leading and subleading corrections in a 1/N S 
expansion of the 2PI effective action: 








(7) 

Figure [T] shows that the NLO quantum correction 
(solid curve) changes the LO result for the distribution 
dramatically. To see how this difference builds up, we 
consider in Figure [2] the total fermion number density 
N^(t) = 8/d 3 pn v ,(t,p)/(2 7 r) 3 as a function of time. 
Both LO (dashed) and NLO (solid) results agree rather 
well at early times at which parametric resonance leads 
to an exponential increase of the inflaton two-point func- 
tion, F$(t, t;po) ~ A^e 2lat . Here 70 denotes the growth 
rate of the maximally amplified momentum mode po 
3]. Significant deviations between the LO and NLO 
fermion production occur at the time where the in- 
flaton field amplitude 4>(i) starts to decrease substan- 
tially, as can be seen from Figure [3l At this time also the 
g<j){t) contribution to the LO equation of motion (|SJ) is 
diminished such that the field independent parts at NLO 
become important. 

Given the somewhat unrealistic parameters for the nu- 
merical results, it is important to discuss what remains 
for realistic scenarios. For this we estimate characteris- 
tic time scales and production rates analytically. Since 
fermions are not yet copiously produced at this time 
can be obtained from the bosonic sector along the lines 
of Ref . Q . For realistic hierarchies this yields 



(8) 



70 \a M 0j 

For physically relevant </>o/Mq ~ 10 6 and erg ~ 0(1) this 



time is indeed of the same order as for the parameters in 
the numerics above. Apart from the fact that tj, depends 
only logarithmically on the mass to amplitude ratio, it is 
well known that for the employed parametric dependence 
of the initial field (po ~ 1/ y/\ the coupling can be scaled 
out of the equation of motion for the rescaled inflaton. 
It remains to estimate the enhanced fermion production 
rate after t^. The amplification at this time is caused by 
the nonperturbativc dependence of on the inverse of 
the coupling. Since the loop correction ((HJ) is ~ g 2 and 
contains a boson line ~ F^ of order 1/A, these together 
lead to corrections of order £ = g 2 /X. From the imag- 
inary part of the corresponding self-energy the fermion 
production rate then scales proportional to £. We have 
explicitly verified this from the numerics for small £ by 
varying g. As a consequence, the shown numerical re- 
sults for g = 0.1 underestimate the even more dramatic 
production for the physically interesting situation with 
£ ~ 1- We emphasize that the observed phenomenon 
of dynamic enhancement of quantum corrections is not 
restricted to chaotic inflation and parametric resonance. 
The large correlations lead to universal behavior [ll[ , and 
the discussion for tachyonic preheating or spinodal insta- 
bilities would follow along the same lines. 

In view of the dramatic differences between LO and 
NLO results, it remains to demonstrate that higher-order 
quantum corrections do not change our understanding. 
Since the interaction term ([1]) is quadratic in the fermion 
fields, we 'integrate out' the fermions on a lattice and 
study the remaining bosonic effective theory by numer- 
ical integration of the equations of motion and Monte 
Carlo sampling [9]. This requires numerical integration 
of the fermion propagator, which becomes feasible in 
3 + 1 dimensions using "low-cost" fermion algorithms 
suggested recently in Ref. [13]. The approach treats the 
fermions exactly. However, the inflaton dynamics be- 
comes classical-statistical, which has been shown in the 
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FIG. 4: Inflaton occupation number of transverse modes. 



FIG. 5: Nonequilibrium fermion spectral function. 



past to be an accurate description for the sufficiently 
early times we are considering here 0, [Tf[ . 

Remarkably, the lattice simulation results shown in 
Figure [T] (symbols) agree well with the NLO approxi- 
mation up to momenta of about 1.3/<7o-Mo [la ] for the 
employed small £. We conclude that the distribution in- 
deed exhibits a thermal value n cq (0) = 1/2 for massless 
fermions at low momenta, which drops at a character- 
istic scale ~ (JqAIq. We emphasize that the inflaton is 
still far from equilibrium at this stage. Figure H shows 
the inflaton occupation number as a function of momen- 
tum at different times, following standard conventions of 
fiefs. [1-3]. In particular, at time t = 250/aoMg the in- 
flaton clearly approaches a nonthermal power-law distri- 
bution ~ l/|p| 4 characteristic for strong turbulence [ill ]. 

A potentially important application of fermion pre- 
heating is the possibility to produce very heavy fermions, 
which may be interesting for the problems of dark mat- 
ter and ultrahigh energy cosmic rays. This question can 
be addressed using the dynamical spectral function ([3]). 
For the discussion of p{x, y) we introduce Wigner coor- 
dinates, X = (x + y)/2 and x — y. For the spatially 
homogeneous system the spectral function then depends 
in addition only on time X . For the temporal com- 
ponent py = tr(7°p)/4 the sum over all frequencies is 
normalized to one, — i J dw/(27r) py(X°; uj, p) = 1, which 
corresponds to the equal-time anticommutation relation 
for fermions. Figure [5] shows NLO results at different 
times and |p|/eroMo — 0. Initially (dash-dotted curve) 
we find a double peak at uj ~ 0, corresponding to mass- 
less fermions, and u> ~ gfio/Nt characterizing very heavy 
fermion states. At later times these peaks become much 
broader and lower, which is characteristic for strong cor- 
relations without a quasiparticle interpretation. 

To conclude, we note that for any boson field with 
nonperturbatively large occupation number and Yukawa- 



type coupling to fermions a similar enhancement of quan- 
tum corrections may be observed. This concerns, in par- 
ticular, also quantum chromodynamics in the presence 
of plasma instabilities. In the context of heavy-ion colli- 
sions for large gluon occupation numbers ~ l/a s , where 
the running gauge coupling a s is small at sufficiently high 
energies, a similar enhancement of quantum corrections 
would have dramatic consequences for quark production. 
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